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Abstract 



Nucleon structure functions, as measured in deep-inelastic lepton scattering, are 
studied within a covariant scalar diquark spectator model. Regarding the nucleon as 
an approximate two-body bound state of a quark and diquark, the Bethe-Salpeter 
equation (BSE) for the bound state vertex function is solved in the ladder approx- 
imation. The valence quark distribution is discussed in terms of the solutions of 
the BSE. 
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1 Introduction 

In recent years many attempts have been made to understand nucleon structure func- 
tions as measured in lepton deep- inelastic scattering (DIS). Although perturbative QCD 
is successful in describing the dependence of structure functions on the squared momen- 
tum transfer, their magnitude is governed by the non-perturbative physics of composite 
particles, and is up to now not calculable directly from QCD. 

A variety of models have been invoked to describe nucleon structure functions. The 
so called "spectator model" is a typical covariant approach amongst them |]TJ. In this 
approach the leading twist, non-singlet quark distributions are calculated from the process 
in which the target nucleon splits into a valence quark, which is scattered by the virtual 
photon, and a spectator system carrying baryon number 2/3. Furthermore the spectrum 
of spectator states is assumed to be saturated through single scalar and vector diquarks. 
Thus, the main ingredient of these models are covariant quark-diquark vertex functions. 

Until now vertex functions have been merely parameterized such that the measured 
quark distributions are reproduced, and no attempts have been made to connect them 
to some dynamical models of the nucleon. In this work we construct the vertex functions 
from a model Lagrangian by solving the Bethe-Salpeter equation (BSE). However, we 
do not aim at a detailed, quantitative description of nucleon structure functions in the 
present work. Rather we outline how to extract quark-diquark vertex functions from 
Euclidean solutions of the BSE. In this context several simplifications are made. We con- 
sider only scalar diquarks as spectators and restrict ourselves to the SU (2) flavor group. 
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The inclusion of vector diquarks and the generalization to SU(3) flavor are relatively 
straightforward extensions and will be left for future work. 

The cross section for DIS of leptons from a nucleon is characterized by the hadronic 
tensor W ,iU (P, q) where P and q are the four- momenta of the target and exchanged virtual 
photon respectively. For unpolarized DIS, the hadronic tensor W^ v is conventionally 
parameterized by two scalar functions F\ and F2. In the Bjorken limit (—q 2 ,P-q^oo; 
but finite x = —q 2 /(2P ■ q)) in which we work throughout, both structure functions 
depend (up to logarithmic corrections) on x only, and are related via the Callan- Gross 
relation: F2 = 2xF\. 

In the scalar diquark spectator model within SU(2) flavor, the valence quark distri- 
butions are extracted from the hadronic tensor (Fig. |l]): 

W^(q,P) = (J + (f)J I Jj^u(P,s)fS(k)Ylm(S(k + q))YS(k) (1) 

xlm(D(P -k))Tu(P,s), 

where the isospin matrix element T3 has to be evaluated in the nucleon isospin space. We 
define u(P,s) as the target nucleon spinor and we use S(k) = l/(m q — ^) and D(k) = 
1 / {mfj — k 2 ) to denote the propagators of the quark and scalar diquark, respectively. The 
integration runs over the quark momentum k, subject to on-mass-shell conditions for the 
diquark and the struck quark. Note that the vertex function T and its PT conjugate f 
consist of two Lorentz scalar functions which depend on k 2 only because the diquark 
is on shell. In the next section we shall determine the vertex functions using a ladder 
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Bethe-Salpeter equation. 



2 Scalar— Diquark Model for Nucleon 

We consider the following model lagrangian: 

C= M0-™ q )A + dr<j>Zd*<p a -m 2 D <l>* a <j> a (2) 
+i^jflMC- l l5 r 2 i, c <f>* a - t -JL=el^ 5 CT 2 ^J </> a , 

where we have explicitly indicated color SU(3) indices only. The symmetric generator 
r 2 of the flavor SU(2) group acts on the iso-doublet field tp for the constituent quark 
carrying an invariant mass m q . The charged scalar field <fi denotes the flavor-singlet scalar 
diquark with invariant mass m^. 

In this model the nucleon with the momentum P and spin s/2 = ±1/2 is described 
by the BS vertex function T (Fig. |^): 

5 ab S(v q P + q)D( VD P - q)tT(q, P)u(P, s) = J p^e* a (0 \T^ a {^)<j> h {--n D x)\ P, s) , 

(3) 

where we set the weight factors to the classical values: r\ q = m q /{m q + m c ) and r] D = 
rri£,/(m q + run)- Then the vertex function V obeys in the ladder approximation the 
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following BSE (Fig.g): 



T(q, P) = g 2 J -Srr^i-k - g - (tfc - Vd)P)S( Vq P + k)D(rj D P - k)T(k, P). (4) 

To solve the BSE for positive energy nucleon states we are free to choose the following 
Dirac matrix structure: 



r(g,P) 



A«(P), (5) 



with A + (P) = 1/2(1 + If/y/P 2 ), the projector onto positive energy nucleon states. 

We assume that the diquark and the nucleon are stable, namely vhd < 2m q and 
\I~P 2 ~ < m q + mo- We then perform the Wick rotation of the relative energy variable 
and choose the nucleon rest frame: P M = P^ = (Po,0). The "Euclidean" functions 
f a (q,P(°>) in terms of the momentum g M = (iq 4 ,q) for a real q A G (—00,00) are then 
functions of qE = V~Q 2 G [0, 00) and z = q^/gs G [—1,1]. We expand each of these 
functions (a = 1, 2) as follows: 

00 

/«(g ) P) = E <B /J?(te) c 'iW. ( 6 ) 

n=0 

where C^(z) are Gegenbauer polynomials and the phase i n is introduced for convenience 
since then the radial functions failE) are real. 
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The BSE in Eq. then reduces to the following system of one-dimensional equations: 



-j 2 oo „ 

A (3=1 m=0 J0 



oo 



dk E IC n S(q E ,k E )W(k E ), 



(7) 



where we have introduced the "eigenvalue" A = (g/47r) 2 for the quark-diquark coupling 
constant. The kernel function )C^(q E , k E ) is a matrix whose elements are real and regular 
functions of q E and k E . By terminating the infinite series in Eq.(^]) at sufficiently high 



state mass, P 2 . To compare the magnitude of the radial functions, let us introduce the 
normalized 0(4) radial functions F n (q E ) and G n (q E ) together with the 0(4) spherical 
spinor harmonics Z n j\ s f2|. In the conventional gamma matrix representation we can 
write the nucleon solution at rest as: 



where q denotes angles for q vector in the four-dimentional polar coordinate system. 
The factor ^/2tt is introduced such that F n (q E ) = fi(q E ), while G n (q E ) is expressed as a 
linear combination of i^ -1 (<?£;) and /^^(Qe)- Thus fi(q, P^) and f2(q, P^) correspond 
to the "upper" and "lower" components of the nucleon Dirac field, respectively. 

Now, let us consider the analytic continuation of f a (q, P^)- Since we are interested in 
applying the BS vertex function to the DIS process, we need to rotate the relative energy 
variable from on the imaginary axis back to on the real one. Recall that the sum over n 



order, we can easily solve Eq. (0) numerically as an "eigenvalue" problem for a fixed bound 




z n F n (q E )Z nh0s (q) \ 
n - l G n {q E )Z nhls {q)l 
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in Eq.(^) converges even for a complex z, if \z\ < 1. We can then analytically continue 
the Gegenbauer polynomials rewriting the argument z — > ( = P ■ qj \/q 2 P 2 for the 
momenta satisfying \(\ < 1. For the radial functions we introduce new functions F n (q E ) = 
F n (q E )/ q E and G n (g|;) = Gn{Q.E)/ 4e 0- We analytically continue these functions by 
changing the argument q 2 E — > — g 2 . We obtain the physical scalar functions f a (q 2 , P ■ q) '■ 

oo 902 /nr /_ n 2\ / \ n-1 

n =i y n(n + 2) J v / 

Note that the Gegenbauer Polynomials together with the square root factors are poly- 
nomials of q 2 , P 2 , and P ■ q, so that each term in the series @ and fllOD is regular and 
real as far as F n (—q 2 ) and G n (—q 2 ) are regular. We may then impose the on-mass-shell 
condition for the diquark and evaluate the sum over n. The resulting f\ and f'2 are then 
functions of the squared quark momentum k 2 and can be applied to the DIS process. 

However, the expressions (|9|) and ( |T0D are valid only for momenta satisfying P 2 \q 2 \ > 
(P ■ q) 2 . Also we found that a naive numerical sum over n, based on Eqs.(||) and flUf ) 
does not converge. Nevertheless, it can be shown that the vertex function T for any 
kinematically allowed k 2 is regular when the diquark is on-mass-shell. This suggests 
that one may be able to continue some appropriate linear combinations of f\ and fi 
outside of this kinematical range. Indeed, we found such a combination which we denote 
by f° n (k 2 ). In terms of this on-shell scalar functions and the quark momentum, k, the 
*The fact that F n and G n are functions of q 2 E was confirmed numerically. 
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vertex function T together with the diquark on-mass-shell condition is given by 



With this on-shell vertex function the valence contribution to the structure function 
Fi(x) can now be calculated from Eq. ([]]). 



In this section we present our numerical results. For simplicity we considered an equal 
mass system; m q = m D = m, and we shall use m as a unit for dimensionful quantities. To 
solve the BSE we used a w-channel form factor. We replaced the quark-diquark coupling 
constant such that g 2 — ► g 2 A. 2 / (A 2 — u) with A = 2 m and u is the usual Mandelstam 
variable. This form factor weakens the short range interaction between the constituents 
and ensures the existence of a discrete bound state spectrum for a large range of P 2 . 

We solved Eq.(|7|) as follows. First we terminated the infinite series in Eq.(^) at 
some fixed value, n ma x- Next we discretized the Euclidean momentum qE and and 
performed the integration over numerically together with some initially assumed radial 
functions /"(fee). This integral generated new radial functions and an "eigenvalue" A 
associated with them. We then used these functions as an input and repeated the above 
procedure until the radial functions and A converged. 

In Fig. £| we plot the normalized 0(4) radial functions, F n (qE) and G n (qE), for the 
bound state mass \f~P 2 ~ = 1.8 m as functions of qs- It is clear that the magnitude of the 




(11) 



3 Numerical Results 
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radial functions with higher 0(4) angular momentum are strongly suppressed compared 
with the lowest ones. This justifies the truncation of the series in Eq.(||). It also confirms 
that the contribution to the "eigenvalue" A from the 0(4) radial functions with n > 4 is 
less than 1 %. The dominance of the lowest 0(4) radial function has been also observed 
in the scalar-scalar ladder model 0. 

In Fig. [| we plot the physical, on-shell scalar functions f™{k 2 ) for the bound state 
mass Vr^ = 1.8m as functions of the squared quark momentum, k 2 . These functions 
are calculated with the maximum 0(4) angular momentum n max = 4. We found that the 
magnitude of f° n (k 2 ) and f% n {k 2 ) are almost the same even for weakly bound states. This 
result suggests that so-called "non-relativistic" approximations, in which one neglects 
the non-leading components of the vertex function (/jj in our model) are valid only for 
extremely weakly bound states: y/P 2 ~ 2m. We have also confirmed that for weakly 
bound states {sfP 2 > 1.8 m) the dependence of f™{k 2 ) on n max is negligible for a small 
spacelike k 2 , e.g., — k 2 < 5 m 2 . However, for large spacelike k 2 , the convergence of the 
sum over n becomes slow for any value of P 2 and numerical results for fixed n max become 
less accurate. 

In Fig.|6| we plot the valence quark distribution F\(x) for a weakly (y/P 2 = 1.8 m, 
solid line) and for a strongly bound state (VP 2 = 1.2m, dashed line). We have used 
n mffl = 4 and the distributions are normalized such that the area below the curve is 
unity. For the weakly bound system, the valence quark distribution has a peak f\ around 
x ~ 1/2. On the other hand, the distribution becomes flat for the strongly bound system. 
TThis peak will shift to x ~ 1/3, if quark and diquark masses such as mo ~ 2m q are used. 
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This behavior however turns out to be mainly of kinematic origin, since the distribution 
function is given by an integral over the squared momentum k 2 carried by the struck 
quark with integration bounds from — oo to fc max , where 



This kinematical bound determines the global shape of Fi(x) to a large extent. 

4 Summary 

We have solved a BSE for a nucleon, described as a bound state of a quark and diquark 
within a covariant quark-scalar-diquark model. We have extracted the physical quark- 
diquark vertex function when the diquark is on-mass-shell from the Euclidean solution. 
The vertex function obtained was applied to a diquark spectator model for DIS, and the 
valence quark contribution to the structure function Fi(x) was calculated. We found 
that the shape of the unpolarized valence quark distribution is mainly determined by 
relativistic kinematics and is independent of the detailed structure of the vertex function. 
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Figure Captions 

Fig. [l|: Diquark spectator process. 

Fig. |2|: The Bethe-Salpeter vertex function. 

Fig. |3|: The Bethe-Salpeter equation in the ladder approximation. 

Fig. f|: The normalized 0(4) radial functions F n and G n for the bound state with the 
mass y/P 2 = 1.8 m as functions of the Euclidean momentum qg. 

Fig. [|: The on-shell scalar functions f° n {k 2 ) (solid) and f'2 n {k 2 ) (dashed) as functions 
of the quark momentum k 2 . The mass of the bound sate is y/P 2 = 1.8m and the 0(4) 
angular momentum is truncated at n max = 4. 

Fig. The valence quark contributions to the structure function Fi(x). The solid 
(dashed) line denotes the weakly (strongly) bound state with the mass: yfP 2 = 1.8 m 
(1.2m). 
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Figure 3: 



14 




Figure 4: 
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Figure 5: 
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Figure 6: 
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